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Abstract 

We calculate the effective Lagrangian for a magnetic field in spinor, scalar and vec- 
tor QED. Connections are then made to SU(Nc) Yang-Mills theory and QCD. The 
Qlh| magnetization and the corresponding effective charge are obtained from the effective 

Q_i! Lagrangian. The renormalized vacuum magnetization will depend on the renormal- 

ization scale chosen. Regardless of this, the effective charge decreasing with the 
magnetic field, as in QCD, corresponds to anti- screening and asymptotic freedom. 
In spinor and scalar QED on the other hand, the effective charge is increasing with 

the magnetic field, corresponding to screening. Including effects due to finite temper- 
ed ' 

ature and density, we comment on the effective charge in a degenerate fermion gas, 
increasing linearly with the chemical potential. Neglecting the tachyonic mode, we 
find that in hot QCD the effective charge is decreasing as the inverse temperature, 
in favor for the formation of a quark-gluon plasma. However, including the real part 
of the contribution from the tachyonic mode, we find instead an effective charge in- 
creasing with the temperature. Including a thermal gluon mass, the effective charge 
in hot QCD is group invariant (unlike in the two cases above), and decreases loga- 
rithmically in accordance to the vacuum renormalization group equation, with the 
temperature as the momentum scale. 
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1 Introduction 



We calculate the one loop effective Lagrangian density for a static uniform magnetic field 
in different gauge theories. From this the effective charge, the magnetization, and other 
relevant objects are obtained. In a recent publication 0, the vacuum of spinor and scalar 
QED was shown to exhibit paramagnetic properties. On the other hand, Nielsen related in 
a well known article ||, asymptotic freedom of QCD to a paramagnetic vacuum of mass- 
less QCD. With the same reasoning, the vacuum of Abelian gauge theories (like QED) 
should exhibit a diamagnetic behavior. In the present letter we resolve the ostensible 
discrepancy between the two different approaches above. This discrepancy originates in 
the use of renormalized quantities in Ref. [|]], verses the use of bare quantities, in order to 
get analogy with a classical dielectricum, in Ref. || . In Section § we shall briefly review the 
properties of a dielectricum, and relate the effective charge to the magnetic behavior. In 
Section | we consider the generic effective Lagrangian for a background magnetic field, and 
how the magnetization is obtained from it. We also give another definition of an effective 
charge, applicable also at finite temperature and density. In Section [| we calculate the 
vacuum effective Lagrangian for spinor, scalar and vector QED in the massive as well as 
the mass-less case. Similar calculations have earlier been performed in for example Ref. || 
(in terms of the vacuum energy), and in Ref. |4[] (in QED). 

Recently much attention has been paid to the /3-function and the effective charge 
in QCD at high temperature ( see e.g. Refs. || [| [F], ||, ||), and high density (see e.g. 
Ref. ]10|]), motivated by the suggested formation of a quark-gluon plasma under such cir- 
cumstances flLTf . The results differ depending on the different choices of gauge, gauge 



fixing parameter J^, and on the vertex considered The effective charge obtained 
from the effective Lagrangian of a background field should be free of these diseases. In the 
background field formalism, as used also in Refs. || [7], f| |9|, the renormalized couplings of 
the different vertices are kept equal and related to the renormalization of the background 
field. We shall in Section |5] calculate the effective Lagrangian in the presence of a heat 
and charge bath, and relate it to the free energy of the plasma. Summing only over the 
physical degrees of freedom, this should be gauge independent. However, to perform the 
explicit calculations we have chosen a certain gauge for the background magnetic field, 
and for the mass-less gluon field. Effective Lagrangians for a static uniform magnetic field 
in a thermal environment have been considered earlier. For scalar and spinor QED in 



for example Refs. [12, 13], and in SU(N C ) Yang-Mills theories (for N c = 2,3) in for ex- 
ample Refs. |TJ, [L5|, [16], [T7|, [18], . However, in the Yang-Mills theories, they focused on 
the effective Lagrangian (or equivalently the thermodynamic pressure) in order to investi- 
gate a possible phase transition, and we disagree with some of the previous results. Since 
perturbation theory in this situation actually only is valid after the phase transition to a 
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quark-gluon plasma has taken place, we shall here mainly focus on the effective charge 
obtained from the effective Lagrangian. A corresponding effective charge has earlier been 
considered in QED |2(| [f3| , but to our knowledge not in non-Abelian gauge theories. 



We shall use a naive real-time formalism, valid here since there are no propagators 
with coinciding momenta, and calculate the thermal effective Lagrangian in spinor, scalar 
and vector QED. SU(N C ) Yang-Mills theory is then related to charged mass-less vector 
bosons, and connections are made with QCD in Section |6|. When possible we compare 
with previous results, and make some corrections. We also consider the effects when a 
thermal gluon mass is taken into account. Finally we discuss the results here obtained, 
and their relevance in Section |7|. 

Since the topics here considered provide an extremely nice and simple example of renor- 
malization, we shall be fairly explicit. 

2 Magnetization and (Anti-) Screening 

In his very pedagogical example, Nielsen described asymptotic freedom (anti-screening) 
in QCD, in the same way as the intuitively clear picture of screening in an ordinary dielectric 
medium. Let us first recapitulate some of the basic features of a classical homogeneous 
dielectricum. The effective coupling of two test charges separated by a distance r is 

where eo is the undressed charge measured in the absence of the medium. The dielectric 
permittivity e(r) must approach unity as r — > 0, since then there is no shielding medium 
between the two test charges. Screening means that the medium will be polarized around 
the test charge, so that the effective charge measured at r > 0, will be smaller than the 
undressed charge (at r = 0). In terms of the dielectric permittivity we thus have e > 1. 
If, on the other hand, the polarization of the medium is such that the effective charge is 
larger than the undressed charge we have anti-screening, and correspondingly e < 1. 

In analogy to this, Nielsen calculated the effective charge in QCD, with the dielectricum 
consisting of the quantum mechanical vacuum, containing virtual particles only. In a 
quantum field theory one will encounter divergences that after a proper regularization 
may be removed by a rescaling of the parameters appearing in the original Lagrangian 
(renormalization). In a classical theory there are no such divergences, so in order to get 
full analogy with the example of a dielectric medium, Nielsen used cut-off regularized bare 
(i.e. before renormalization) quantities. An ultraviolet momentum cut-off A, corresponds 
to a smallest distance tq = 1/A, where the bare charge is measured, corresponding to the 
classical undressed charge, at r = 0. 
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The relativistic invariance of the vacuum of a quantum field theory requires that the 
permittivity is connected to the magnetic permeability /i pcrm , through 




(2.2) 



This has no counterpart in an ordinary polarizable medium. It turns out that it is easier 
to calculate the magnetic susceptibility x (such that fi pcrm = 1 + x) i* 1 a background 
(color) magnetic field B, than to directly calculate the permittivity in a (color) electric 
field. By a heuristic reasoning, Nielsen related the distance to the field strength, according 
to r « 1/yeB (we are assuming eB > 0). We thus find the dielectric permittivity 



In terms of bare regularized quantities (denoted by the subscript "6"), we thus find that 
screening (e& > 1), corresponds to Xb < 0, i.e. diamagnetism; and anti-screening (ej < 1), 
corresponds to Xb > 0, i.e. paramagnetism. 

How could then possibly the Abelian gauge theories discussed in Ref. exhibit a 
paramagnetic behavior, when they cannot be asymptotically free? 

The answer to this lies in the renormalization procedure. We can never measure such 
things as bare quantities, but renormalization is necessary in order to obtain the physical 
parameters. The renormalization is performed at some momentum scale A, arbitrary but 
finite. The reference charge will be the renormalized charge, measured at momentum 
scale A, or equivalently on a distance n ~ 1/A. For the clarity of the reasoning, and 
the comparison with the classical dielectricum, we are assuming some sort of on-shell 
renormalization, where the renormalized charge corresponds to the charge measured at 
the renormalization scale. In general it is sufficient only to remove the divergences in the 
renormalization procedure. Then one will have some finite relation between the physical 
charge actually measured, and the renormalized charge that just is a parameter of the 
theory. In analogy to the classical dielectricum we thus have 



where e\ is the charge measured at distance n, i.e. £\{t\) = 1, corresponding to the 
charge renormalized at momentum Ai ~ 1/Vi. The general criteria for screening is that the 
effective charge should decrease with the distance or increase with the momentum scale, 
and vice versa for anti-screening. When considering distances smaller than r 1; screening 
instead corresponds to E\ < 1, and anti-screening to E\ > 1. Without knowledge of the 
scale at which the reference charge is measured, the magnitude of e, will thus not tell 
whether we have screening or anti-screening. 




(2.3) 




(2.4) 
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In the renormalized spinor and scalar QED considered in Ref. [|1] the renormalization 
is performed at the lowest momentum scale (as is standard), corresponding to r\ — > oo. In 
this case we thus always have r < n, and screening here corresponds to e < 1, which means 
paramagnetism x > 0, and similarly anti-screening would correspond to e > 1, and thereby 
X < 0, i.e. diamagnetism. However, if the theory exhibits anti-screening (like QCD), the 
effective charge is becoming infinitely large at large distances, so the renormalization cannot 
be performed at vanishing momentum scale. 



3 Effective Lagrangian and Effective Coupling 



The generic effective Lagrangian for a background magnetic field may be separated into 
its different contributions 

•^eff ^tree £v&c ^mat • (3-1) 

Here 

£tree = —B 2 /2 , (3.2) 

is the free tree level part; £ vac is the vacuum contribution due to virtual particles; and £ mat 
is the contribution due to real particles at finite temperature and density. The vacuum 
contribution £ vac is calculated to one-loop order in Section £|, and the matter contribution 
£ ma t will be considered in Section ^. 

From the effective Lagrangian we may obtain the magnetization 

d 

M EE M vac + M mat = ^(£vac + £mat) • (3.3) 

The vacuum magnetization, originally proposed in Ref. is a real physical quantity, 
but only gives measurable effects at extremely high field strengths. Performing another 
derivative we find the magnetic susceptibility 

d 

X = Xvac + Xmat = ^ M w + M mat ) . (3.4) 

on 

At vanishing temperature and density we may use the vacuum magnetic susceptibility to 
obtain the effective charge according to Eq. ( |2.3| ) and Eq. fl2.4| ). In the presence of matter 
at finite temperature (T) and chemical potential (/i), Lorentz invariance is broken when 
choosing a preferred frame of reference in which the medium is at rest and in equilibrium. 
Then there is no connection between the permittivity and the permeability. We must 
therefore find another way of obtaining the effective charge in this case. It may be done 



through the identification [20 



£ eff - ~\bI, = , (3-5) 
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where we have used that eB = e e gB e g is required to be invariant under renormalization. 
Performing the derivative with respect to (eB) 2 , that is invariant under renormalization 
as well as Lorentz transformations (B 2 — E 2 is Lorentz invariant and here reduced to B 2 ), 
we find 

e cS ~ d(eB) 2 e 2 eB e ' { } 

where M is the magnetization. Notice that one may also define the magnetic susceptibility 
as the response function of the magnetization M = \B. In this case we find 



111 1 1 M 



e 2 e 2 1 + x e 2 eB e 
to the lowest order in the coupling. 



(3.7) 



An external magnetic field H is introduced by adding a term £ ext = j^ xt A u to the 
effective Lagrangian. Neglecting a surface term we find £ ex t = B • H. By construction 
£ e fj is invariant under renormalization. If C c g + C cxt is to be invariant, the invariance of 
eB requires H/e to be invariant. This seems quite reasonable, since V x H = j ext . The 
external field H is thus only a function of the external charges, and thus should scale as a 
charge. Minimizing £ e ff with respect to B we find the mean-field equation 

B = H + M vac (B) +M mat (5) . (3.8) 

This equation is telling us how to find the average microscopic field B, that also is the 
acting field felt by the particles in the medium, in the presence of an external field H. 



4 The Vacuum of Spinor, Scalar and Vector QED 

We shall here consider the one-loop vacuum contribution to the effective Lagrangian of 
spinor, scalar and vector QED with an external static uniform magnetic field B = Be z . In 
order to ease the physical interpretation, and the comparison with Ref. ||, we shall use a 
cut-off regularization procedure throughout this Section. 

The Lagrangian for a spin 1/2 particle of charge — e coupled to the electro-magnetic 
field is 

£ ferm = -~F^F^ + ^-m)V> , (4.1) 

where the covariant derivative is D u = d v — ieA u , and we have used the shorthand notation 
p= 7 M -D^. The corresponding Lagrangian for the spin case reads 

£SC ai = _ p^F^ + Dl<jfD u <p - m 2 ^(f) . (4.2) 
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For reasons that will become obvious in section |6] we are also interested in the case of a 
spin 1 vector particle with gyro-magnetic ratio 7 = 2. The corresponding Lagrangian is 
then obtained by adding an anomalous magnetic moment term to the minimally coupled 
Lagrangian |[21||, with the result 



£ vec = ~F^F» V -^(D^Wl- DlWl)(D»W u - D u W»)+m 2 WlW» (4.3) 

-ze(7 - — e . 

The general quantized theory described by C vec is not even renormalizable [p!q| . But here 
we shall consider no virtual photons. The renormalizability then follows in exactly the 
same way as for scalar and spinor QED. We separate the vector-potential into A u + A u , 
where A u is the vector-potential for B, and A u corresponds to the quantized radiation field. 
We shall here only calculate the one-loop effective Lagrangian. We may then neglect the 
radiation field, apart from the kinetic term —F 2 /2. For simplicity we shall not write out 
the radiation field, but only add the contribution from thermal photons (or gluons) in the 
end. 

A particle of charge — e, spin cr and gyro-magnetic ratio 7, has the spin magnetic 
moment fi = — / ~fe(cr) /2m. With spin projection s along — eB the energy spectrum reads |3[] 



E n .s = ^m 2 +p 2 z + eB(2n + 1 - 7 s) , (4.4) 

where the dependence on the momentum in the direction of the field (p z ) has been sup- 
pressed, and we assume eB > 0. 

It is illustrated how to calculate the vacuum contribution at one-loop level in spinor 
QED in Refs. 0, |l3 |, and also in scalar QED in Ref. [|13|]. In the general case we find 



similarly, using the trick to perform the derivative with respect to the mass 

^ = j( _ ir+1 ^ ? £j^_^_ , (4 . 5) 



s n=0 ' 



where a is the spin of the intermediate particles. We have here suppressed that the contribu- 
tion at B = should be subtracted, due to the normalization of the generating functional. 
The subscript b denotes bare quantities, i.e. before renormalization. The difference in sign 
comes from the fermionic Grassmann algebra. 

We may now use Cauchy's theorem to perform the integration over u. Integrating with 
respect to m 2 we find, with the subtraction explicit 

eB 00 r 

£vac,b = (-1) 2 ' +1 7^E E J d P* E n,s ~ £vac,b(5 = 0) . (4.6) 



s n=0 ' 



In order to interpret this expression we shall calculate the density of states for particles 
in a box V = L x L y L z . We then need that the wave functions in each of the theories here 



6 



considered, with the choice of gauge A 

1 



(0, 0, Bx, 0), are of the generic form 



n,py,p z ,s' 



x,t) 



2^/2E„ 



exp[±i(-E njS t+p y y+p z z)} I n 



I'n 



x 



7T / 



cxp 



--eB (x-2l 
2 V e5 



(4.7) 



(4.8) 



where if n is a Hermite polynomial. Without the magnetic field, the number of states on 
the momentum interval dpi is Li/2n dpi. With the above magnetic field, this still holds true 
for i = z. The wave-functions in Eq. (fO§) corresponds to a harmonic oscillator, centered 
at Xo = Py/eB, that has to be within the box, i.e. < p y /eB < L x . Since the energy is 
independent of p y we must sum over all possible values. The resulting degeneracy of states 
is then 



V 



eB 



(4.9) 



The vacuum Lagrangian density is thus equal to the negative vacuum energy density 

-1 



vac,b 



V 



[E va , c — E va _ c (B 



(4.10) 



again without the contribution at vanishing magnetic field. The vacuum energy was the 
starting point in Ref. [§]. The vacuum energy needs to be regularized. As it stands it is 
quadratically divergent, but when the B = part has been subtracted the result is only 
logarithmically divergent. 

Instead of immediately integrating with respect to m 2 , we shall here first use 

jE = ^Q^- E ^ dX + ^ ' < 4 ' U > 

where we have introduced the ultra-violet cut-off A, that essentially removes the contribu- 
tions for E > A due to the exponential suppression. We may now perform the Gaussian 
integral over p z , sum the infinite geometrical series in n, and integrate with respect to m 2 . 
Subtracting the contribution for B = 0, and changing variable of integration to t — x 2 , we 
arrive at 

eBt 



(-l) 2a „ /-oo dt . 9 \ 

£vac,b = : R 2 2^ / ,Trexp(-mt)< . 



exp(eBt"ys) — 1 



(4.12) 



For eBi^/s — 1) > m 2 this is divergent for large t. This is the case for large fields in vector 
QED (er = 1) with 7 = 2, that will be treated separately below. Let us now first consider 



the massive case in Section 14.11, and then the mass-less case in Section 



4.1 The Massive Case 

In the limit as x = eBt — > 0, we have 



x 



— 7— exp(7sa;) = 1 + ^sx + \x 2 (~f 2 s 2 - - 
smh(x) 2 V 3 



0{x 3 ) 



(4.13) 
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Notice that X) s s = 0. Define from the quadratic term 



X = 



16tt 2 



? (- v " I) 



(4.14) 



Let us now subtract and add the term containing 



/~ a f exp(- m 2 t)=EU-VA 2 ) = -7, + ln(-^)+o(^) , (4.15) 

where E\ is an exponential integral, and 7^ = 0.57721566 ... is Euler's constant. In order 
to get rid of the cutoff A we must now perform a renormalization. The Ward identities 
of the theories here considered requires the product eB to be invariant. Let us therefore 
rescale the external field and the coupling according to 



B b = Z]! 2 B{\) , 



(4.16) 
(4.17) 



where e(A) and B(X) are the charge and background field renormalized at momentum scale 
A, respectively. Adding = ln(A 2 /A 2 ), we may write 



C 



cS 



'2^-2 (ZA - 1) ^ + 2 (e5)x 



-7£ + In i 



'A 2 / 
A 2 



+ C vac (eB,X) , (4.18) 



where the finite, renormalized vacuum contribution is 



C vac (eB,X) = heB) 2 X \n( ^) + £ v UeB) , 



C wac (eB) 



2' ' " \m 

^ exp( - mt Mw 



-1 



,2(7 



E 



exp(e_Bt7s) — 



eBt 



sinh(efit) 



- 1 



- x- 



(4.19) 
{eBt) 2 1 



We must now choose Z\ in such a way that the divergence as A — > 00 is removed 



Z x -1 = e 2 (A 2 )x 



-7b + In ( ~\2 



A 2, 
A 2 



(4.20) 



(4.21) 



i.e. 



e 2 (A) 



1 

12 +X 



Performing the derivative with respect to A we find the lowest order /5-function 

de(X) 



X- 



dX 



(3[e(X)\ = -xe'(X) 



(4.22) 



(4.23) 
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In the case of (spinor) QED, a — 1/2 and 7 = 2, we find 



X ferm = ^ • (4.24) 

Similarly in scalar QED, a = gives 

X scal = K • (4.25) 

48vr 2 v ; 

Inserting this into Eq. ( |4.23| ), we recognize the correct /3-function of spinor and scalar QED. 
In the theory of vector QED with gyro- magnetic ratio 7 = 2, the /3-function really does 
not exist, due to the lack of renormalizability. However, for vector bosons interacting only 
with the external field we find 

7 

v vec | = — A 26) 

A I massive \Q7l^ ' ' 

This does not agree with the renormalization of the coupling considered in Ref. |22[ (corre- 
sponding to x vcc — 5/167T 2 ). Notice the crucial difference in sign of \ for vector QED, that 
we later shall relate to QCD. The solution to the above renormalization group equation 
may be written 

With x < this corresponds to the effective charge increasing with the momentum scale 
A, whereas x > corresponds to the effective charge decreasing. Obviously e(A) and B(X) 
are defined such that £ e ff is independent of A. We have 

£eff = + \{eBf X In(^) + LUeB) . (4.28) 

Notice that Z\ only depends on A through the dimension-less parameter A/A. As long as 
the A dependence is absorbed in Z\ (as is necessary to get finite expressions) the RGE fl4723|) 
will not depend on how Z\ and thus e(A) are defined. Since the cutoff A no longer appears 
in the effective Lagrangian, we may now send it back to infinity where it belongs. 

In order to obtain a physically clear picture let us now consider the high field limit. 
Substitute y = eBt in Eq. (^4.20|) , and split the integral at y , such that eB/m 2 3> y ^> 1. 
The contribution from y < yo is easily seen to be 0[(eB) 2 ]. The contribution from y > yo 
is dominated by the subtracted term x(e£>t) 2 /2, that gives 

£vae 1 , m 2 r°° dy ( m 2 \ 1 / eB\ 

~ ~nX / — eX P EV --nX h K • ( 4 ' 29 ) 



(eB) 2 2^ (eB) 2 J yo y * \ eB a ) 2^ \m 2 

At the energy scale eB = A 2 ^> m 2 , the two logarithms will cancel, and thus the effective 
Lagrangian assumes a purely Maxwellian form 

C cS ~ ~B 2 (\ 2 = eB) , e£ > m 2 . (4.30) 
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We thus find that at least in the limit A>m, that e(A) and B(X) are the charge and field 
strength measured at momentum scale eB = A 2 , respectively. 

In the high field limit (eB > m 2 ), Eq. ( ^29|) gives to leading order 

M V ac(A) . (eB\ 

Xvac(A) = -e 2 (A)xln(^) . (4.32) 

We see that the sign of the magnetic susceptibility in addition to x depends on the relative 
magnitude between eB and A 2 . This does anyhow correspond to screening in spinor and 
scalar QED (x < 0); but anti-screening in vector QED (x > 0) since 

dXvac{X) - -e 2 (A)x4 ■ (4-33) 



d(eB) y e£ 

The permittivity e is thus decreasing with the energy scale eB in spinor and scalar QED 
(with x < 0) , but increasing in vector QED (with x > 0). 

In the high field limit we find the effective charge 

' 1 +xln{ ^) , eB^>m 2 , (4.34) 



e 2 eS (eB) e 2 (A ) * V A 2 

that is a solution to the lowest order RGE with the momentum scale identified as A = y/eB. 
Notice that e e s must be independent of A by its definition in Eq. (|3.6|), since C e s and eB 
so is, and that this follows from Eq. ([4.27|) . It is no coincidence that the above effective 
coupling satisfies the lowest order RGE, since it corresponds to a summation of the leading 
logarithms. The high field limit is dominated by the term containing x, subtracted in the 
renormalization, that is defining the /3-function. 

The effective Lagrangian will assume a particularly simple form if we choose A = m 

£ eff = ~B 2 + C vac (eB) . (4.35) 



Here we have suppressed the dependence on A, since in the weak field limit 

(4.36) 



1 

2' 



C cS = -1B 2 + B 2 



' c 2 (eB) 2 
m 4 



so that e and B are the ordinary electric charge and magnetic field, measured at vanishingly 
small magnetic field, corresponding to long wavelengths. The vacuum of spinor and scalar 
QED will thus be paramagnetic (Xvac > 0) in the high field limit {eB 3> A 2 = m 2 }. 
Actually, in the spinor case it is paramagnetic for all values of the magnetic field, as shown 
in Ref. [|IJ, to which we refer for a more extensive treatment of the vacuum magnetization 
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in spinor and scalar QED. In spinor QED we have the vacuum contribution to the effective 
Lagrangian explicitly 



eBt coth(eBt) - 1 - \{eBtf 



In scalar QED we find 



seal _ 
vac " 16 



1 r°° dt . 2 , 



siuM^)- 1 + 6 (ei?t) 



(4.37) 



(4.38) 



As remarked above the case of vector QED for eB > m 2 needs more careful considera- 
tions. Let us use the similarity between vector and scalar QED to write 



fvec 
•^"vac,b 



massive 



3£^ b + A£^ b . (4.39) 



Written on the vacuum energy form corresponding to Eq. (f4.6|) the extra term reads 

eB 

Obviously, this will contain an imaginary part for eB > m 2 



A ^e,b = -7^a / d Pz (^m*+pl-eB-ie - y/m*+p* + eB ) . (4.40) 



P R r VeB-m 2 , (eB—m 2 )eB 

1* AC™ = — / d Pz ^eB-m 2 - P 2 = Q(eB- m 2 ) { - — J— . (4.41) 

Working out the real part in analogy to the general case we find 
l{eB)\ ( A 2 \ eBm 2 




2 27T 2 \^m 2 + eB)\m 2 -eB\) ^ 

-C{eBf + iQ(eB - m 2 ) ^ ~^ eB , (4.42) 

where C is an irrelevant constant that anyhow will be renormalized away. Notice that the 
result also may be obtained by analytical continuation from the result valid for m? > eB, 
i.e. by substituting \m 2 — eB\ \— > (m 2 — eB — is). Also notice that the dependence on the 
cut-off A is independent of the magnitude of m 2 — eB, so the renormalization will not be 
affected. The imaginary part in the effective Lagrangian is telling us that the configuration 
of a spin 1, gyro- magnetic ratio 2 particle in a strong (eB > m 2 ) static uniform magnetic 
field is unstable. This problem has been addressed earlier in the literature, see Section ^| 
for references and a discussion. The conclusion is that a condensate will be formed in the 
unstable mode. The magnetization of this condensate will alter the background magnetic 
field so that it is no longer uniform on a microscopic scale. The energy of the lowest mode 
will then become positive so that the instability is removed. However, the net change in the 
magnetic field is found to be small, so we simply neglect this here. We are thus assuming 
our results to be valid also in vector QED, and neglect the appearance of the imaginary 
part. The coefficient \ in front of [eB) 2 ln(A 2 /m 2 ) obtained from Eq. (|4.39 ) agrees with 
Eq. CT . 
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4.2 The Mass-less Case 



Let us now consider the mass-less limit of spinor, scalar and vector QED. For vector bosons 
it is then necessary to choose the gauge D^W^ = 0, in order for the wave-functions to be 
of the same form as in Eq. (|4.7|) . This will actually introduce Faddeev-Popov ghosts. In 



for example Ref. |T7| it was shown that the contribution from these ghosts exactly cancel 
the contribution from the unphysical degrees of freedom. This means that if we restrict 
the sum over polarizations to s = ±1, we do not need to consider ghosts. 

In this mass-less case we may not subtract the next term in the expansion of Eq. (|4.13|) , 
since this would cause a divergence for large t. Instead we substitute x = eBt in Eq. (|4.12|) 
and integrate by parts to find 



C 



^, = \x(eBflJ^) + C -(eBf 



(4.43) 



where the u.v. finite constant is 



<,7T Z 



roo d 

dx In x — 

, Jo dx 



exp(7sx) 
x sinh x 



(4.44) 



In the mass-less limit, vector QED will be unstable for all values of B. Since only two 
different polarizations s = ±1 now are possible, we write 



C 



vec 
vac 



9£scal , A £ 



vec 
vac 



(4.45) 



We may directly take the mass-less limit in Eq. Q4.42 ). The result is that Eq. ft4.43|) still 
is valid, but with an imaginary part in the constant C . This is of no surprise since the 
renormalization is due to ultra-violet divergences, i.e. p 2 z ^> eB. Again we shall neglect 
this imaginary part and the instability. The effective Lagrangian of Eq. ( |4.43|) is what was 
considered in Ref. (but in terms of the vacuum energy) . Nielsen extracted from here 
the leading bare vacuum magnetic susceptibility 



Xvac,b ^ e^xln 



A! 

eB 



For mass-less vector bosons we find explicitly 
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24^2 



(4.46) 



(4.47) 



This has the opposite sign as compared to spinor and scalar QED. Bare vector QED (that 
we in Section |^ will relate to QCD) will thus be paramagnetic, that corresponds to anti- 
screening and asymptotic freedom. However, performing the renormalization we have 



-- x (eB) 2 ln — 



(4.48) 



12 



The sign of the renormalized vacuum magnetic susceptibility, i.e. if we have a paramagnetic 
(Xvac > 0) or diamagnetic (Xvac < 0) vacuum, will thus again depend on the relative 
magnitude between the field strength eB, and the renormalization scale A 2 . 

Asymptotic freedom really means anti-screening with the effective charge vanishing at 
vanishing distance (infinitely large momentum scale). From the RGE this is obvious since 
e(A) = is a fixed point. We may also compare the effective charge at two different scales 



-eff 



[eB) 



-cff 



(eB 



eB 

l + 4f(e£o)xln( — 



(4.49) 



If x > 0, as in vector QED and QCD this quotient is vanishing as eB — > oo, i.e. we have 
asymptotic freedom. In spinor and scalar QED with x < 0, the coupling appears to grow 
infinitely large at a finite (but extraordinary large) value of the magnetic field. However, 
perturbation theory (we are only considering one loop effects here) cannot be extrapolated 
that far. 

5 Thermal Spinor, Scalar and Vector QED 

We shall in this section consider the contributions to the effective Lagrangian due to finite 
temperature and density of particles. We may on this one loop level use the naive real 
time formalism. In Ref. it was shown that the substitution in the fermion vacuum 
propagator 

/fM= 2 L^. ^-2n5(u 2 -E 2 )f F (u) , (5.1) 
uj z — tj A + IE 

in order to obtain the thermal part of the propagator, works also in a magnetic field. Here 
the fermion one-particle distribution in thermal equilibrium is the Fermi-Dirac distribution 

where fi is the chemical potential, and 1/(3 = T is the temperature. In the bosonic case we 
substitute similarly 

2 L, ■ ^27r5(LU 2 -E 2 )f B (uj) , (5.3) 
where Jb, in thermal equilibrium, is the Bose-Einstein distribution 

, , . e(w) e(-w) 



Performing these substitutions in Eq. (4.5), we may use the 5-function to integrate over u. 
Integrating with respect to m 2 , the constant of integration is determined by the fact that 
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the exponential suppression from the particle distributions requires /3mat — * 0, as tti — > oo. 
We then find the contribution from the heat and charge bath to the effective Lagrangian 



eB 00 r 
£mat= (2^??Si dp 



s n=Q ' 



Pi 



1 



1 



e P{E n , a -ti) _ ( — 1) 2ct e P(E n ,s+v) — ( — 1) 2 ° 



(5.5) 



Integrating by parts with respect to p z we recognize the free-energy density 



£mat = lri Z 



:-i) 



2cr+l 



eB 



(3 (2tt) 2 



00 . 

EE/^ {in [l - (-l) a *e-«*"-'0 



s n=0 



+ ln[l-(-l) 2CT e- /3( ^^ t) ]} , 
where Z is the partition function of the gas. Let us now split £ ma t into 

•^mat = •^-'mat,0 -£mat,l j 

where £ ma t,o is the field independent part 

£mat,0 



(5.6) 



(5.7) 



1 )2T+1 £ /'-^{liM (_i)^ e ^^(p)-M) 



/3 



(27T) 



In 



2a-l3{E{p)+ l i) 



(5.8) 

whith -E(p) = + P 2 - We are particularly interested in the high temperature expansion 
of the effective Lagrangian. We shall consider the massive case for spinor, scalar and vector 



QED in Section |5.1| , and investigate the corresponding mass-less limits in Section |5.2| . In 
every case the contribution from thermal photons 



9tt 2 

,-phot _ fJ_rnA 
^mat — ^ 



(5.9) 



should be added to the effective Lagrangian. In Section ^73] we consider also a degenerate 
fermion gas at low temperature and high density. 



5.1 The Massive case 



The case of massive fermions in a magnetic field was treated in Ref. ||13|| . Here we just 
quote the final result 



/'ferm 
■'-mat,! 



duj®(uj 2 — m 2 )f F (uj) 



\ 1 f°° tt_ -tf?P-m 2 ) 



\4vr 5 /2 7 fB/2 



[eBtcoth(eBt) - 1] 



1 j,/< U\ :i 2 . fir 



u 2 — m 2 ' 



, / (-T-) sin I irl 

2vr 3 ^V I J V 4 eB 



(5.10) 



The high temperature, weak field limit {T 2 ^> m 2 ^> eB; \l = 0} was considered in 
Ref. HI!, with the result 



CaT + 4r = ^ - ix fcrm (e5) 2 ln (^) + 0[(eB) 2 } 



(5.11) 
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Notice here that in the high-temperature effective Lagrangian, terms with higher powers 
of eB are suppressed as (eB) 2 (eB/T 2 ) n ~ 2 , n > 4. To leading order there are thus no 
non-linear electro-magnetic interactions in high temperature QED, in agreement with a 
diagrammatic analysis |p3|. Also notice that the lnm 2 terms will cancel when this is added 
to the effective Lagrangian in Eq. ( |4.28| ), and that the effective charge obtained from this 
will satisfy the lowest order RGE with the scale A = T. 



Let us now consider the case of massive scalar QED. For the sake of completeness we 
here quote the result (see e.g. Ref. [ ]E|] ) 



/"•seal 
mat,l 



JduQ(u 2 -m 2 -eB)f B (u)l^ J - 2 J o 



An 3 



5x3/2 



sin 



7T 

4 



7rZ 



u) 2 — m 2 



00 dt 
-eB" 



j—t(aj 2 —m' 2 ) 



eB 



eBt 



smh(eBt) 



(5.12) 



In the high temperature, weak field limit, this gives 

r 2 

45~ 2 
T (eB) 2 



/^scal I /'seal 



^T 4 - U(eB) 2 



In 



m 



m 487T 



l + O 



IE 



o 



+ o 



m f eB\ 
T 



(5.13) 



In the massive vector case we again write L 
we have 

dPzP 2 z 



vec 
mat 



Af vec 

^''-mat 



eB 



1 



q /'seal I A fvee 

"'-'mat T l-ii-mat 



1 



£ -i e^ -- 1 - 1 E 00 e 



/3£ ,o 



T 2 J 

For m 2 > eB and fi = 0, 
' (5.11) 



The calculations performed in order to find the high temperature expansion of this expres- 
sion are explicitly shown in Appendix |A.1|. The final result reads 



A rvee 
z - xz -'mat 



(eB) 2 
4tt 2 

eBm 2 
' 8tt 2 



In 



In 



T 2 (4tt) 



m 4 - (eB) 2 \ 
m 2 + eB \ 



eBT - 

+ —— (Vm 2 + eB 
2n 



eB\) 



eB 



2tt 2 



2 + eB\ 3/2 



+iQ(eB - m 2 ) 



eBT^feB 



2tt 



eB[eB-m 2 ) 
8^ 



(5.15) 



Notice that the ln(m 2 ±eS) terms cancel between Eq. ( |4.42| ) and Eq. ( |5.15| ). For eB > m 2 
the lowest energy mode will become unstable, resulting in the imaginary part as depicted 
above. The result follows also in this case by analytical continuation m 2 — > m 2 — is in 
expression valid for eB < m 2 . Notice that the T-independent imaginary part exactly 
cancels the imaginary part of the corresponding vacuum contribution in Eq. fl4.42|) , and 
that the resulting imaginary part is increasing linearly with the temperature. 
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Of course such an imaginary part is not acceptable. Since a full treatment of this 
imaginary part is out of the scope here we are left with two alternatives. Either neglecting 
the imaginary part, or neglecting the total contribution from the tachyonic mode, that also 
has a real part. In the former approach we may immediately use Eq. ( |5.15 ). The latter 



approach is considered in the mass-less case in the following Section |5.2| . The massive case 
(that is relevant for spontaneously broken gauge theories) can be treated in full analogy. 

5.2 The Mass-less case 

The explicit calculations performed in order to obtain the high temperature expansion 



in mass-less spinor and scalar QED are presented in Appendix |A.2|. Adding the thermal 



contribution for B = 0, £ ma t,0j The final result reads in the spinor case 

+ 47a? = 4 - \x lem {eBf In ^ + 0[(eB) 2 ] . (5.16) 

Just like in the massive case, all non-linear electro-magnetic interactions from the vacuum 
part has been cancelled by contributions from the thermal part. In the scalar case we find 

+ Cat = ^ - ^^IC(-^)l(e5) 3/2 T - In ^ + 0[(eB) 2 } . (5.17) 



Again the ln[(ei?) 2 /A 2 ] has been cancelled by thermal contributions, but there is in this 
case also the term linear in T. 

We are now left with the vector bosons. The configuration of mass-less vector bosons in 
a static uniform magnetic field is unstable for all field-strengths. We will here neglect the 
contribution for p 2 < eB — m 2 in the lowest Landau level in order to avoid the imaginary 
part. However, this may affect the high temperature behavior, that is not solely governed by 
large momenta, but also by small energies when the Bose-Einstein distribution is becoming 



very large. The final result of the calculations performed in Appendix |A.3| reads 



AC 



eB) 2 [ T 



mat o 1 / 

7T VeB 



IT 2 n 
- In + 1 

2 4eB 2 



2 



+ \ ln 7B +0{1) \ ■ (5 ' 18) 



The same expression is obtained if we substitute m 2 — > — is in Eq. (|5.15|) , and subtract the 
corresponding integral f^** analytically continued to imaginary energy. 

5.3 The Dense Degenerate fermion Gas 



Chodos et al. |24] has suggested that the fermion matter contribution in the limit of large 
chemical potential may be of importance in heavy- ion collisions. Let for simplicity T = 0. 
The leading behavior of the first term in Eq. ( |5.10|) is then 

ix fcrm (e5) 2 ln^ 2 
2 m / 



£ reg ^--X fcrm ( eJ B) 2 ln^ . (5.19) 
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We may write the oscillating second term in Eq. ( |5.10| ) as 
r , 

*— 'nsr 



(eBfl 2 ™ 1 t \ 

— -L^ Im ex P 



Here we have defined 

Jn = 



duj exp —iirn—- 



( ix m 2 \ 

i nn — 

I 4 eB ) 



(5.20) 




—mn 


u 2 ' 
eB_ 








jerf 




— erf 


■ u nnm 2 
e W4 / 


c 


V e£ 


V eB 



(5.21) 



where error functions were identified. Let us now neglect the second term that is indepen- 
dent of the chemical potential. Using the asymptotic expansion of the error function Q 
for 7r/i 2 ^> eB we find the leading term at large chemical potentials 



(eBf' 2 g 



4vr 4 /i 



n=l 



5/2 



Sill 



n 



it a 

h nn 

4 eB 



m 



(5.22) 



This is suppressed at large \i, but performing the derivative with respect to eB we find the 
leading behavior, with mod[v4] = A — int[v4], 



/i 



e 2 B 



1 00 1 

—Y 



n=l 
1 



3/2 



sin 



2n7r 



;U 2 — m 2 ' 



2e5 



rC 



•-,mod 



/i 2 — m 2 
2eB 



(5.23) 
(5.24) 



In Ref. 01 t ne results of Ref. [|13|] was used to find Eq. ( p.23| ), using another method. 



Chodos et al. 0] suggest that the effective charge obtained from this, increasing linearly 
with n for fi 2 = m 2 + 2eBn, n G Z, could be of relevance in heavy ion collisions. In 
Ref. 0] it was claimed that this is valid only for fi 2 ^> eB ^> m 2 , but here we have shown 
that this is the leading behavior for fi 2 ^> eB, m 2 irrespective of the relative magnitude of 
eB and m 2 . The effective charge then seems to become divergent as eB — > 0. This must 
be an artifact due to one out of two possibilities. 

1. The break-down of perturbation theory, as the coupling is becoming stronger. 

2. The derivative in the definition of the effective charge. If we instead would define the 
effective charge by l/e 2 ff = — 2C e s/(eB) 2 , then the contribution from the oscillating 
part is suppressed at large chemical potentials. 

Furthermore, the linearly increasing amplitude in M osc is a result of the sharp de Haas-van 
Alphen oscillations at the Fermi surface, that will be smoothed out at finite temperature. 
The (-function in Eq. (|5~23D takes its minimal value ((-1/2,0) = £(-1/2,1) ~ -0.208, 
and its maximal value ((—1/2, 0.3027) ~ 0.0934. The effective charge is thus oscillating, 
and we doubt its physical significance in this limit. 
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6 QCD with a Background Magnetic Field 



Quantum Chromodynamics with a background magnetic field has been extensively studied 
in the literature. It is well-known that the mean-field equation (|3.8|) has a non-vanishing 
solution for B, even in the absence of an external field (i.e. H = 0) in QCD. The vacuum 
energy is lower in this state than for B = 0, as pointed out in Ref. f25| . However a 
tachyonic mode appears 
potential 



cf. Eq. (O), that causes an imaginary part in the effective 
signaling the instability of this configuration. This unstable mode has been 
suggested to be removed by a (1+1) dimensional dynamical Higgs mechanism [26||. The 
corresponding condensate shows a domain-like structure, and may form a quantum liquid 
of magnetic flux tubes |27], 28 1. A corresponding treatment has also been done in the 
electro- weak theory (e.g. Refs. |29|, |30|). Here the change in the magnetic field due to 



the condensate, that assumes a lattice structure, is found to be small compared to the 
uniform field. It has also been suggested in another approach that the tachyonic mode 
could be stabilized by radiative corrections |31, p2[ , or by the condensation of an auxiliary 
. In Refs. 133, 



field 



the imaginary part was found to appear also for a non-Abelian 



like background field. However, it was argued |34] that the imaginary part originates in 
the abuse of the formula C e s oc In det G _1 , where G is a propagator, valid only for positive 



definite G. It was concluded 34 that the contribution from the unstable modes and the 



imaginary part should not be trusted. In Ref. |19| the imaginary part was shown to vanish 
for large enough values of a color condensate A Q . 

Since a full treatment of the tachyonic mode is out of the scope of this monograph, some 



approximations or assumptions are required. In for example Refs. |L4|, [T5| the contribution 
from the tachyonic mode was neglected, as we have done here in Section |5.2| . However, in 
Refs. [|i~6| , |i~7H it was argued that the contribution of the unstable mode should be taken into 
account by analytical continuation, as in Eq. ( p,15| ). Starting from the partition function or 
the free energy of the quark gluon plasma, cf. Eq. ( |5.6D , it seems very unnatural to include 
an unphysical tachyonic mode. Due to the presence of the thermal distribution function 
the thermal contribution of the tachyonic mode will not be purely imaginary, unlike the 
vacuum contribution. In Ref. |36| the existence of a zero-energy mode, after the removal 
of the tachyonic mode, was pointed out. This suggests that including contributions only 
from p 2 z + m 2 — eB > in the lowest energy mode is a reasonable approximation. However, 
this will exclude the contribution from the plausible condensate. We shall therefore briefly 
comment also on the corresponding results when the real part of the contribution from the 
unstable mode is included. These problems may be solved by introducing a thermal gluon 
mass, large enough to remove the instability. The effects of such a mass are considered in 
Section |(T4| . 

Usually a temperature dependent (as well as RGE running) coupling is obtained from 
the vertex correction. In the background field formalism here employed the vertex cor- 
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rection is related to the vacuum polarization. We shall therefore consider the vacuum 
polarization, and review the obtained effective couplings in Section O . 

The conserved charge in QCD may correspond to the color charge, the baryon num- 
ber or the quark flavor. Each quark is carrying the color charge "1" of some color, and 
the baryonic number 1/3 that relates the baryon chemical potential to the color chemical 
potential. Nature only seems to allow for equal amounts of the different colors, that corre- 
sponds to equal chemical potentials. However, the gluons carry equal amount of color and 
anti-color (but not the respective), resulting in a vanishing chemical potential for them, as 
well as their linear combination in terms of the W fields. The W bosons carry an equal 
amount of color and anti-color, in the ideal combination for interactions with the external 
field. The problem with Bose-Einstein condensation in the lowest energy mode thus never 
occurs, since \x = for the gluons and W bosons. Nevertheless, we shall here mainly focus 
on the high temperature situation with vanishing charge density (i.e. fi = 0). Since a finite 
chemical potential only affects the quarks, we may in this case immediately use the results 
obtained for QED in Section |5.3|. 



6.1 Connections with Spinor and Vector QED 

We shall here first relate the theory of mass-less spin a — 1, gyro- magnetic ratio 7 = 2 
bosons interacting with a background magnetic field, to SU(N C ) Yang-Mills theory in a 
background chromo-magnetic field (cf. Refs. 0, P7[). The Lagrangian of SU(N C ) Yang- 
Mills, is 

£yk = -7 E G%G^ = E £ {a) , (6-1) 

4 0=1 o=l 

where 

Gl v = d,El-d v El + gY,F bc ElEl . (6.2) 

b,c 

We shall now consider this theory in the presence of a background chromo-magnetic field, 
that we choose in the a' = N% — 1 direction in color space, i.e. 

Ef'-^A. + Ef^ . (6.3) 

Here is the vector potential corresponding to a static uniform magnetic field, and E^ -1 
is the quantum field. As in the Abelian theories considered before, we are interested in 
the one-loop effective Lagrangian for this background magnetic field A^. To the one- loop 
order we may neglect every occurrence of g, when not in the combination gA^. We then 
find 

1 1 7V c 2 -2AT?-2 

= —F^ - -gF^ Y, E r' 3k E^Et , (6.4) 
4 1 j=i k=i 
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where F^ v = d^A u — d u A^ is the (color) electro-magnetic field-strength tensor. As in the 
Abelian theories considered above we have not explicitly written out the radiation field 
E^~ l , but we must remember to include its thermal contribution. In the limit g = 
the Lagrangian of Eq. ( |6.1| ) corresponds to (N^ — 1) "photons", that gives the thermal 
contribution for B = to the effective Lagrangian. For f ^ a' we find 



£<f) 



iV?-l 



d.El - d v E? - g £ r' j ' k (A,E* - A u E k ) 



k=i 



(6.5) 



a'j'k 



where we have used the total anti-symmetry of the structure constants p ' a k = — f 
Generally we may choose the generators of SU(N C ), such that f a 'i' h is non- vanishing only 
for one k = k', for fixed a' and f. Let us now define 



K - ^[^ i ' fc,) +^°"' fc ' )t ] 



E 



In terms of the charged vector boson field W we have 



where we have defined the covariant derivative 

Adding the relevant term from C^ a '^ in Eq. (|6.4|) , we find 



(6.6) 



+ C (k') = _i[ D * w y,w _ D * w u',k>)^ D » w u',k>)v _ D v W w,k>)^ ; (67) 



(6.8) 



£(f'.*0 --F 2 = -If 2 - - \D*W {j '> k ' )] - D*W ij '^]\D fl W {j '' k ' ) ' y - D u W ij '' k ' )fl] 

^ a 2 ^ A 4 v v u ] L - 



W ' fj, 

—igf F^ ; - - - 



(6.9) 



Comparing with Eq. (|4.3|) , we see that this is exactly the Lagrangian for a spin a = 
1, gyro-magnetic ratio 7 = 2 boson of charge — gf a ^' k , interacting with a background 
field, described by A^. The terms in the original Lagrangian £ym describing the self- 
interaction of the E a fields, that we here have discarded when considering the one-loop 
effective Lagrangian for A u , are essential for the renormalizability of the theory [j37fl . The 
total Lagrangian £yii * s then obtained by summing over all such pairs (j', k') 



C 



No 
YM 



--F 2 



ihk) 



(6.10) 
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where (j, k) are the pairs with non-vanishing f a '^ k , for a' = N% — 1. Again using the total 
anti-symmetry of the structure constants, the total charge squared of these particles is 



JV?-liV?-l 



2 IV j . 

L = s 2 E(f' J1 ) 2 4E E/ a '^r' jfe • (6.11; 



(i,fe) 2 i=i fc=i 

Now the SU(N C ) generators in the adjoint representation are (T® ec )jk = f a ^ k . We thus find 
a group invariant squared charge, independent of the direction a' in SU (N c ) 

elJ ah = yTr [T^J = jN c 6 ab . (6.12) 

We now wish to couple fermions in the fundamental representation to our Yang-Mills 
theory. The corresponding Lagrangian reads 

£& = W + jEI*J?-^ , (6-13) 

0=1 

where the color indices j = 1,2, . . . , N c have been suppressed. In the background field 
of Eq. ( p.3|) , we again neglect g when not combined in gA^, and find 

£ i ™ = m + gAT f i m -m)^ . (6.14) 

Generally we may choose T f g' rm diagonal, i.e. T^ Tm = diag(ti, ti, . . . , tjvj- Comparing with 
Eq. ( |4.1|) , we see that the Lagrangian in Eq. ( |6.14|) corresponds to iV c fermions with charges 
e = —g(t 1 ,t 2 , . . . , ijVc)> coupled to the external field described by A^. The squared sum of 
their charges is then also group invariant 

eLJ ah = L 2 E ^ 6ab = 9 2 ^ PSoAJ = 9 ^ ab ■ (6-i5) 

Invariance under background field gauge transformations requires the product eB to be 
invariant under renormalization (see e.g. Ref. 0). We may then immediately use the 
results of the previous Sections. We shall for simplicity use the relations between quadratic 
charges in Eqs. (|6.11| , |6~1~5|) . When terms not quadratic in the coupling appear in the 



effective Lagrangian, the correct approach is to sum over the moduli of the constituent 
charges to the power considered, as we assume eB > 0. However, when not only containing 
even powers of the coupling e, the result is not group-invariant. This means that it depends 
on the specific direction in color space chosen for the magnetic field. Any such result is 
dubious, and probably unphysical. 
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6.2 The Effective Coupling in Thermal QCD 



In the previous Section we related QCD with a background (chromo-) magnetic field to 
the corresponding cases in spinor and vector QED. However, for the contributions to the 
effective Lagrangian independent of the magnetic field, we must instead identify and sum 
over all degrees of freedom. We shall here consider Nf flavors of quarks with possibly 
different masses m/. They all come in N c different colors. The number of gluons is N 2 — 1. 
Since we only are summing over the true degrees of freedom, either in the vacuum energy, 
or in the free energy, no Faddeev-Popov ghosts are needed. The effective Lagrangian for 
SU(N C ) Yang-Mills with a background color magnetic field, and Nf flavors of quarks in 
the fundamental representation may thus be written as 

,2 



r QCD 
*"eff 



£ tree +iV c Y:£^ + (iV c 2 -l)/: 



vec 
mat,0 



'ferm I 2 
\ dr. 



/ , ^"matj^ferm 
/ 



mat,l V vec 



(6.16) 



with m i— > rrif in the different terms in the sum over flavors. In vacuum the corresponding 
renormalization gives the correct lowest order QCD /^-function 



A^ = /3 QCD [s(A)] 



X QC V(A) 



(6.17) 



where we obtain from Eqs. ( 4.24 , 4.47 ) and Eqs. ( 6.12 . 6.15| ) 
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(6.18) 



In the high temperature limit {T 2 ^> eB, m 2 } for all flavors of quarks, we find when 
neglecting the contribution from the tachyonic mode 
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(6.19) 



Here the ordinary renormalization at T = 0, /i = has been performed at the renor- 
malization scale A. The dependence of A in g (A) is cancelling the explicit A dependence, 
so that the effective Lagrangian is independent of A, as follows from its definition. No- 
tice that the {gB) 2 ln(m 2 ) terms cancel between Eq. (|4.28 ) and Eq. ( |5.11| ), so that there 



is no dependence of the different quark masses whatsoever, as long as T 2 3> m 2 , and 



that the same result follows from the mass- less case in Eq. ( |5.16|) . Due to the appear- 
ance of (gB) 3 / 2 it is not quite correct to use the squared average charge for the glu- 
ons. In SU(3) with B in the iV 2 — 1 = 8 direction in color space, we have the relevant 
structure constants / 458 = Vo/2 = f 67S (see e.g. Ref. |H]]). We should thus substitute 
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(iV c /2) 3 / 4 = (3/2) 3 / 4 ^ 2 (3/4) 3 / 4 in the term containing (gB) 3 / 2 T, but this depends on the 



direction in color space. In Refs. JTj, fL5|l , where the tachyonic mode was neglected, no high 
temperature expansion was made, and the graphs presented do not extend to temperatures 
high enough, that we could compare our result with theirs. 



Let us now use Eq. ( |3.6|) to define a field and temperature dependent effective coupling 
g e s(T,eB). Comparing the effective charge at different temperature and field strengths, 
we find 



1 



1 



&(T,eB) g 2 eS (T ,eB ) 
where we have defined 

F(x) =x[y 



+ x QCD In -2 + F(T/eB) - F(T/eB ) 
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(6.20) 
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This effective coupling is thus decreasing as a function of the temperature. The leading 
behavior, with x = T / \fe~B is 



9eA X ) 
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(6.22) 



This is a faster decreasing coupling than predicted by the vacuum RGE with the identifi- 
cation A i— ► T. The dominant T/y/eB ln(T 2 /eS) term is neither present when considering 



the running coupling obtained from the vacuum polarization (see Section |6.3j ), nor when 
including the real part of the contribution from the tachyonic mode as considered next. In 
this latter case we find instead 



pQCD , ^QCD 
•'—mat,! ' •'—vac 



2tt 



/Q\ 3/4 

(J [(l+z) + 2(^2-l)C(-l/2)] 
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(6.23) 

We have here summed over two vector bosons with the charge e = g^/3/2, to obtain the 
coefficient in front of (gB) 3 ^ 2 T for comparisons with earlier results. The factor 2(3/4) 3 / 4 
should not be present in SU(2) Yang-Mills theory (with / 123 = 1). 

Evaluating numerically we find [1 + 2(y / 2 — 1)£(— 1/2)] ~ 0.827781, that exactly equals 
the corresponding numerical coefficient in Ref. |~7| (where SU(2) was considered). As a 
matter of fact, using the same identification with a generalized Riemann ( function as used 
here, the coefficient of Ref. [0 may be written as i+ [2 3/2 |C(- 1/2, 3/2) | -1] ~ i+0.827781. 
However, we do not agree with the high temperature limit in Ref. |16|| , and particularly 
not with the sign of the imaginary part. That sign of the imaginary part is obtained if we 
substitute m 2 — eB ie — eB (i.e. opposite to the Feynman prescription in the propagator) 
of Eq. (|5.15|) . It was pointed out already in Ref. [0 that the different sign on the imaginary 
part is unphysical, and corresponds to a blow up instead of a decay of the corresponding 
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configuration. In Ref. |19) (where SU(2) was considered) a coefficient of (eB) 3 ^ 2 T/2n in 
perfect agreement with Eq. ( |6.23| ) is obtained (in the limit Aq = 0, that is what here is 
considered), but the real contribution from A£^ t (i.e. "1" in 1 +i) is unfortunately lost 
in the final result. Probably a typographic error also has caused a coefficient in front of 
the {gB) 2 ln(T 2 /A 2 ) not in accordance to the RGE in Ref. Considering the real part, 
Eq. (|6.23|) gives an effective charge that is increasing with the temperature. The leading 
behavior of F in Eq. (|6.20|) is in this case 

/S\ 3 / 4 3 r- 
F{x)=z2(jJ -[l + 2( v / 2-l)C(-l/2)] • (6.24) 



Notice that the effective coupling in Eq. ( p. 22 ) and the one obtained from Eq. ( 6.24 ) are 



not invariant under transformations in color space. In Eq. ( |6.24|) we have explicitly stated 



the result when the magnetic field is chosen in the "8" direction in SU(3) in order to 
compare with previous results. We feel dubious to the physical significance of such an 
effective charge. 

6.3 The Vacuum Polarization 

At finite temperature, the broken Lorentz invariance results in two different possible tensor 
structures in the vacuum polarization 

Il^(k) = P% u Tl T (k) + P£ u Tl L (k) , (6.25) 

where P^ v and P£ v are the spatially transverse and longitudinal polarization operators, 
respectively. Considering the effective coupling obtained from the vacuum polarization for 
a gluon with momentum scale |k| = k in mass-less thermal QCD, we may write 

1 1 ■ ^QCD ln ^ + fi(T/K) _ n(T /K ) • (6.26) 



g 2 (n,T) g 2 (K ,T ) wjj 

We have here already used the well known form of the vacuum polarization in the absence 
of matter, leading to Eq. (|4.27|). Chaichian and Hayashi || suggest to use either g 2 Il = 



KJ* 2 = {k 2 / k 2 )ICJ k 2 , or g 2 fl = IIIJk 2 = n^'at - (k 2 /K 2 )UZ t )/(2K 2 ). We shall 
here only consider the static momentum configuration, k° = 0, |k| = k. 

With the formalism used in this work we may easily calculate n™ at in the limit k v — 0. 
Using the similarity between the way the chemical potential /x and the vector potential Aq 
enters the Lagrangian, we find 

n^° at (^ = 0)=e 2 ^^ . (6.27) 
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Since we have found the quark masses to be irrelevant in the high temperature limit, let 
us for simplicity consider the mass- less case only. Performing the derivative in Eq. ( |A.9| ), 
we may then let /i = 0. The field independent part is then obtained as 



n°° at0 = ^-T 2 YY v v 9 /""^exp[--l . (6.28) 

mat.O 27r 2 Z-Z^ p J Q t 2 ^ 2V V ' 

This gives = n^ato™ = e 2 T 2 /3, and in the vector case n^J|o = 2e 2 T 2 /3. What 
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roodt . 1 



amounts to the field dependent part we proceed much in the same way as in Appendix |A. 2 . 



Substituting u = (3 2 l 2 eBt/2, we may directly perform the Poisson resummation, since the 
I = term is vanishing in this case. Again it is necessary to separate out the k = term 
for bosons. In the scalar case this gives a contribution 



.OU.sci.l '^7Tm2 — V^2 r 



n^:r = -Ve5T^-[-C(l/2)] . (6.29) 



The other terms are found to be suppressed at high temperatures. Identifying the momen- 

tL _ n oo 

mat mat 



turn scale k 2 = eB, we thus find in QCD, using for k = II mat = n°° 



n(x) = x 2 (^ + ^)-x^^?|C(l/2)| + 0(l/x) . (6.30) 

The leading T 2 behavior agrees with Ref. JJ, but the coefficient in front of x = T/k, 
approximately 0.1365A ;r c /2, does not. The absence of InT in Eq. ( |S.30|) indicates that the 
identification between the magnetic field eB and the momentum scale k 2 does not work 
here. 

On the other hand, using the gauge invariant Vilkovisky-De Witt effective action in- 
dicates that it is the transverse part of the vacuum polarization tensor that governs the 
renormalization ||. In this case a cancelation occurs between the leading T 2 /k 2 terms. 
The coefficient in front of the next leading T/k term is found to depend on the gauge fixing 
parameter. However, the Vilkovisky-De Witt effective action speaks in favor of the Lan- 
dau gauge £ = 0. Recently Elmfors and Kobes used a Braaten-Pisarski resummation 
scheme to calculate II^ at self-consistently in the high-temperature limit. Their result 
reads in our notation 

n(x)^-x^£[(3 + 2 + 14] . (6.31) 

This will give an effective charge increasing with the temperature. The effects of the 
inclusion of a non-perturbative magnetic mass was also considered in Ref. |§. In the 
Landau gauge this did not qualitatively change the asymptotic behavior of the effective 
charge, whereas it could do so in other gauges with £ > 1. Recently, Sasaki |39| has 



used the pinch-technique [|(| to calculate a gauge- invariant thermal /5-function in QCD. 
It was explicitly shown in various gauges that the thermal /3-function is invariant, and 
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the result agrees with the one obtained in the background field formalism in the Feynman 
gauge £ = 1. No Braaten-Pisarski resummation was performed in this case, but we have 
found no reason why Eq. ( |6.31|) with £ = 1 should not be the resummed pinch technique 
result. We have no answer why the Vilkovisky-De Witt effective action, on general grounds 
supposed to be gauge-invariant, gives a different result from the pinch technique, explicitly 
shown to give equal results in various gauges. 



6.4 Inclusion of the Thermal Gluon Mass 

It is well-known that the solution of some infra-red problems in thermal field theories 
requires a resummation of the dominant diagrams to all orders in perturbation theory |38 



The result is concluded in the elegant "Hard Thermal Loop Effective Action". Here we 
are primarily interested in regularizing the infra-red behavior. We shall therefore only take 
the thermal gluon mass into account. The derivation performed here will only be heuristic. 
Considering the high temperature limit, we shall use the well-known result |38| for B = 



m 2 G = u 1 (k 2 = o)^ [ n c + ^-)9y ' (6 - 32) 

that must be the leading contribution also for T 2 ^> gB. This is the on-shell self-energy 
for transverse gluons only, corresponding to the polarizations s = ±1 for the charged W 
bosons. We thus add and subtract a mass term m 2 G /2 J2a ^a-^u t° the SU (N c ) Lagrangian 
in Eq. ( |6.1| ). The subtracted mass term is necessary in order not to change the original 
Lagrangian, that would ruin gauge invariance. It is treated as a counter term, in order that 
the hard thermal loops are not counted twice. The resummed transverse gluon self-energy 
is for example n^ s = n T — m G . The contributions from this counter term are calculated in 
Appendix [FJ. Rewritten in terms of W, the resummed Lagrangian corresponds to Eq. (|4.3|) 
with m 2 = m 2 G . Therefore, we may immediately use the results from massive vector QED, 
but only with two polarizations (s = ±1) similar to the mass-less case, since a dynamical 
mass term cannot change the number of degrees of freedom. Furthermore, we shall assume 
the thermal mass to be so large that the instability is removed, i.e. m G > gB. With these 
assumptions we may immediately use our previous results in the high temperature limit 
to find 



2^ w ' A 2 m G 2 48vr 




(6.33) 



where g and B denote the charge and field renormalized at the momentum scale A, re- 
spectively. Actually we should here also have included the thermal quark mass m 2 = 



26 



[(N 2 — l)/2N c ]g 2 T 2 , but this gives no relevant effects. Notice that expanding according to 
gB < m G in the last term of Eq. (|6.33|) will only give even powers of 
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(6.34) 



G 



so that the result is group invariant to any order. The correct approach would be to use 
e vec = ^ 4 /2 Tr [T" e ' c T" e ' c T" e ' c T" e ' c ] etc. for higher powers. Here we shall only consider the 
leading behavior for tuq ^> gB. The contribution from the counter terms is then 
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(6.35) 



We may now perform the derivative with respect to B to find the magnetization 
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We have here written "12 — 1" in the last term to indicate that this coefficient enters exactly 
as x vcc = 2% scal + A% vec . In terms of the charge renormalized at momentum scale Ao, we 
now find the effective temperature dependent charge 
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(6.37) 



Notice that mc oc T, so that the only T dependence is in the logarithm. The other 
constant terms may be renormalized away. The effective charge obtained with a thermal 
mass regularizing the instability is thus decreasing logarithmically in accordance to the 
vacuum RGE with A = T. This was the original assumption by Collins and Perry |jTT|, on 
the formation of a quark-gluon plasma (a notion later coined) as a result of asymptotic 
freedom. Also notice that the result follows provided that the thermal mass is mp oc T, 
and m/i > gB, regardless of the explicit form of mp. We could therefore have used the 
electric mass m 2 E = H L (ko = 0, k) = 3wig, or the (nonperturbative) magnetic mass m 2 M = 
n T (^o = 0, k) = (cg 2 T) 2 . The thermal mass used here follows also in the uniform limit 



U T (k ,k = 0) 



m 



g- 



We believe that it is the thermal gluon mass used here that is 



relevant for the physical polarizations of the gluons. 



7 Discussion 



In this paper we have calculated the effective charge of different gauge theories in vacuum 
and in a thermal environment. In the vacuum case the effective charge is related to the 
magnetic susceptibility. The general criteria for (anti-) screening, is that the effective 
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charge is (increasing) decreasing with the distance or the inverse momentum scale. If 
the reference charge that we are comparing our effective charge with is taken as the bare 
charge, measured at infinite energy, the effective charge is always smaller than the bare 
charge in spinor and scalar QED. This corresponds to screening, and is related to bare 
diamagnetism. In QCD, on the other hand, the effective charge is always greater than the 
bare charge. This implies anti-screening, asymptotic freedom and bare paramagnetism. 
However, in spinor and scalar QED the reference charge is customary taken as the classical 
charge, measured in the infinite wave-length limit. In this case the effective charge is 
always greater than the reference charge, and the corresponding vacuum magnetization 
will show a paramagnetic behavior. These theories do anyhow exhibit screening since the 
effective charge is increasing with the energy scale. In QCD there is no such natural scale 
at which to define the effective charge, since the charge is becoming infinitely large in the 
long wavelength limit. The magnitude of the effective charge compared to the reference 
charge , and thereby also the sign of the magnetic susceptibility, will therefore depend on 
the relative magnitude between the scales at which the reference charge and the effective 
charge are measured. 

However, some caution is required in the interpretation of asymptotic freedom in QCD 
in terms of anti-screening in a dielectricum. In the presence of quarks and gluons, external 
fields are screened also in QCD, see e.g. Ref. |^TJ. Chromo-electric fields are screened 
with the electric (Debye) mass m% ~ g 2 [(N c + N f /2)T 2 /3 + £/ /xj/(27r 2 )]. To cure some 
infrared problems also magnetic fields are believed to be screened with a non-perturbatively 
generated magnetic mass m 2 M = 0(g 4 T 2 ). We believe that this could be related to the 



condensate removing the tachyonic mode, as suggested already by Cornwall |)40| , but this 
needs to be investigated further. We have here found that the anti-screening is caused by 
the large spin magnetic moment of the gluons that themselves carry the color charge, in 
terms of the W field in Eq. (|6.6| ). We may view the anti-screening as an effect of dispersion 
of the color charge jlT]. For example, a static blue (B) quark may become red (R) by 



emitting a BR gluon. This will effectively distribute the blue charge over a volume ~ r 3 . 
When investigating the charge in a volume A~ 3 <C r 3 , only a small fraction of the net blue 
charge is found. 

In the thermal case we mainly focus on the effective charge in QCD. There are several 
advantages in this approach of calculating the effective charge from the effective Lagrangian 
in a back-ground magnetic field. Background field gauge-invariance is maintained, even 
though we only consider one particular choice of gauge here. This gauge-invariance implies 
that the product of the coupling and the background field are invariant under renormaliza- 
tion. The quark-gluon, three gluon and four gluon couplings are thus renormalized in the 
same way, and kept equal. Summing over physical degrees of freedom only, in the vacuum 
energy and the free energy, there is no need to introduce ghosts. The running coupling 
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obtained from the RGE is a function of scaled momenta — *> e 4 /^, but this means that we 
move off-shell. Most often the effective charge is calculated in the static limit fc M = (0, nk), 
relevant for the screening of external fields. This means that one is in the deep Euclidean 
region k 2 = —k 2 . The running coupling so obtained cannot be ascribed to real particles, 
only to virtual particles in internal processes. The effective coupling obtained from the 
effective Lagrangian, on the other hand, is directly related to the interactions of the real 
particles in the heat and charge-bath. Moreover, the thermal contribution to the effective 
Lagrangian is related to the free energy. Since physical, measurable quantities are obtained 
by performing derivatives on the free energy, it must be gauge invariant and cannot depend 
on the gauge fixing, up to an irrelevant constant. 



At fixed magnetic field, we find from Eq. ( |6.20| ) the leading behavior of the effective 



coupling in QCD as a function of the temperature 
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As found when explicitly calculating parts of the vacuum polarization tensor, we cannot 
transform the dependence on the magnetic field, to dependence of momentum through 
k, 2 rs eB. However, we may consider the limit of vanishing magnetic field. In the limit 
gB — > in Eq. ( [7.1| ) we find the simple behavior 



9lM) T 



(7.2) 



i.e. a coupling linearly decreasing with the temperature. The effective charge obtained in 
the limit of vanishing magnetic field has the advantage of being independent of the direction 
in color space for the magnetic field, unlike the effective charge in Eq. (|7.1|) as indicated 
by the non-analytical behavior in terms of N c . Notice, however, that the appearance of 
the (g_B) 3 / 2 Tln(T/A) term in the effective Lagrangian, that enforces asymptotic freedom 
at high temperatures, is solely due to the negligence of the tachyonic mode. If this mode is 
taken into account, by analytically continue m 2 — > —ie in Eq. ( |5.15| ), this term disappears. 
But we are in that case left with an imaginary part in the free energy. Considering the 
real part only we find in the high temperature limit an increasing effective charge 

n 2 (T m 9ls(To,gB) 
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Also this effective charge suffers from non-invariance under transformations in the color 
group SU(N C ), and thus cannot be trusted. In this case we cannot take the limit gB — ► 0, 
since this would take us across the Landau pole, at which the coupling is becoming infinitely 
strong. We obtain qualitatively different results whether the tachyonic mode is neglected, or 



29 



the real part of its contribution is taken into account. However, particularly when starting 
from the free energy, we consider it physically more reasonable to neglect all contributions 
from the imaginary energy mode, that of course cannot be present in a real situation. The 
group invariant charge obtained in the field free limit when neglecting the contribution 
from the tachyonic mode could thus be taken as an indication of the true situation. 

A thorough investigation of this tachyonic mode has so far required numerical treat- 
ment, and the (small) change in the magnetic field due to the condensate in this mode 
would also alter the other modes present, and thereby also the effective Lagrangian. This 
seems to destroy the advantage of simplicity in obtaining the effective charge, using the 
effective Lagrangian in a magnetic field. However, including a thermal gluon mass may 
take us out of this dilemma. We find in the high temperature limit when the thermal 
mass is large enough to completely remove the instability, an effective charge decreasing 
according to the vacuum RGE (|4.27|) with A t— > T 

1 1 / T x 2 



This is what naively has been expected for long |Tl[ , but hitherto not obtainable. Obviously, 
this effective charge is also valid in the absence of a (chromo-) magnetic field. However, 
the derivation outlined here is merely heuristic, and more rigorous considerations are to be 
be performed. 

If we instead consider the case of large flavor chemical potentials, we may immediately 
use the expansions in Eq. (|5.19| ) and Eq. fl5.23| ). It is the latter, oscillating contribution 



that is most interesting. It is oscillating with an amplitude \M osc \/e 2 B oc n/yeB, but it 
may be positive as well as negative depending on the ratio (/i 2 — m 2 )/{2eB). The effective 
coupling thus may increase as well as decrease. However, due to these rapid oscillations we 
are doubtful to this definition of the effective charge in the limit of large chemical potentials 
and low temperatures. 
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A High Temperature Expansions 



A.l Massive Vector QED 

We shall here find the high temperature limit of 
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we may perform the Gaussian p z -integral in Eq. ( |5.14| ). We can also identify ^2 
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where K v = K_ v is a modified Bessel function. We may now write 
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Using the derivative property of Bessel functions 43 , we find [42 



This finally gives 
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In order for this to be valid also for eB > m 2 , we must consider the analytical continuation 
defined by m 2 — > m 2 — ie. The same result is obtained also if we treat the contribution for 
pi < eB separately in analogy to the vacuum case. The result is presented in Eq. ( p,15| ). 
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A. 2 Mass-less Spinor and Scalar QED 

In the mass-less case, we have found it necessary to add the vacuum contribution in or- 
der to perform a Poisson resummation and find the high temperature behavior. We have 
not managed to perform this analysis using the cut-off regularization previously utilized. 
Therefore, we shall in this section instead use a dimensional regularization in 4 — 28 di- 
mensions. In order to keep the coupling dimension-less, we substitute e t— > \ s e, where A is 
an energy scale. For E UjS >/iwe may use the expansion in Eq. ( |A.2j ), and the identity in 
Eq. ( |A.3| ). Performing the summation over n, we find 
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Integrating over p z , subtracting the B = part £ mat ,o and substituting u = j3 2 l 2 eB\ 5 t/2, 
we find 
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We now wish to perform a Poisson resummation, and must therefore add the I = term 
to extend the summation to Ya^-oo- This term with I = is the bare vacuum Lagrangian, 
regularized in 4 — 25 dimensions, cf. Eq. ( |4.12 ). Let us now for simplicity consider the 



case of a neutral plasma, /i = 0. Again the vector case needs special care, so here we first 
consider only spinor and scalar QED with a = 1/2 and a = 0, respectively. Performing 
the Poisson resummation we find 



C vac + C M = (eB)-^=^_j EJ Q ^exp 
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The term with k = for bosons (a = 0) gives a contribution to C^ t , finite for 5 = 

{eBf/ 2 T r du I u \ {eBfl 2 T /2 1 

^rV^io ^jsinhR" 1 / = ^ 2_( 2/ " 4 )v^C(-2) , (A.12) 

where £(—1/2) ~ —0.207886 , and we have used T(— 1/2) = —2-^77. We have here identified 
Riemann's zeta function ((z) from an integral representation in Ref. E| that is valid for 
Re (z) > 1. However the subtracted "1" on the left hand side of Eq. (|A.12|) is regularizing 
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this expression (the integral is manifestly convergent), and we have numerically checked 
the equality. For the other terms in the sum, we sum over s and expand [f|2] 

1 = o^-^L 7^1 B *nu , (A.13) 



8vr 2 [ sinhfw] J 2 A 8tt 2 ^ (2n)\ 

for 7cr = 0, 1. We may then integrate over u. Due to the alternating sign of the Bernoulli 
numbers B 2n , the modulus of the sum is smaller than the modulus of the first neglected 
term. For v > 2 we may let 5 = 0, and find the contribution 

1 2 2n - 2(1-7(7) , N9 / eB \ 2( ^ x) r(2i/-3/2) >/ s , A s 

that is suppressed for large T. From the first term in the expansion of Eq. (|A.13|) we obtain 
a contribution 

where Yl'k means that the term with k = should be excluded for a = 0. We can now 
identify a generalized Riemann's (-function [ j^j 

1 00 1 

E'^pg ^ 2g [k + (1 _ g)]1+M ^ 2Cd + M, 1 - a) . (A.16) 



As 5 — > we have the expansion [[12 

((1 + 25,1-^ = 1+7 + 2(7^2 + 0(5) , (A.17) 

for a = 0, 1/2. We now expand r(§+5) = ^+0(5), {X 2 /ttT 2 ) s = l-51n(7rT 2 /A 2 ) + C(5), 
perform a renormalization to absorb the divergent term 1/5 and all terms oc {eB) 2 in the 
bare coupling, and then let 5 vanish. 



A. 3 Mass-less Vector QED 

In this case we cannot employ a technique similar to Eq. ( |A.12j ), since the integral would 
become divergent due to the exp(7sw). Here we shall neglect the tachyonic mode and start 
integrating at p 2 = eB — m 2 in the contribution from the lowest Landau level. Since the 
lowest energy then is vanishing, only /i = is possible in this case. If we have a finite charge 
density it must therefore reside in a Bose-Einstein condensate in the zero energy mode. 
When the tachyonic mode is not included, integrations by part will cause non-vanishing 
surface terms. We shall here start from the free energy in Eq. (|5.6| ), that we believe is most 
physical. We then have 

A£™ = - / d ft ln[l - e~^\ - \ d V M\\ - e-^°.°] . (A.18) 

7T A UVeB JO J 
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Substitute x = (3E 0jS , and split the integral at 1 ^> xq yeB/T. Expanding for x > xo 

+ . (A.19) 



x 2 + eB/T 2 Jx 2 -eB/T 2 T 2 x 3 \T 2 J 



The remaining integral is convergent, so this will only produce an irrelevant 0(eB) 2 term, 
and terms suppressed at large T. For x < xq we instead expand 

lnfl-e"*] =\nx-^ + 0{x 2 ) . (A.20) 



We may now perform the integrals over x to find Eq. ( |5.18 



B Thermal Gluon Mass 



We shall here briefly outline the results from the inclusion of a thermal mass m@. Let us 
start without the matter contribution. Considering one-loop effects only, we may integrate 
out the particles in the generating functional of Green's functions. Adding and subtracting 
the term corresponding to the thermal mass, we find in the simplest case of a mass-less 
(i.e. without the thermal mass) scalar field 

J d 4 xC cS = J d 4 x (-^ 2 ) + i Tr ln[i(-D 2 - + mj)] . (B.l) 

We now wish to define our perturbative expansion in terms of the massive field. Let us 
therefore split the logarithm and the effective Lagrangian according to 



J d 4 x (£ vac + 4i) = i Tr ln^-D 2 - m% )] + i Tr In 



1 + 



D 2 — ml 



(B.2) 



Performing the mi derivative on the first term, we recognize the propagator (— D 2 — m|) _1 . 
This may then be generalized to arbitrary spin, resulting in Eq. ( |4.5| ) for the case of a 
magnetic field considered here. Expanding the second logarithm in Eq. ( B.2 ) to leading 



order, we again find the trace of the propagator. Similar generalizations to arbitrary spin 
give 

p R 00 r 1 



Integrating over u, we may then use Eq. ( |4.11|) . In this counter term we shall encounter 
no divergences, so we may immediately let 1/A = 0. We may then perform the Gaussian 
integral in p z , and sum the infinite geometrical series in n. Substituting t = x 2 , and 
subtracting the contribution for B = 0, we find 
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Expanding for eBt < eB/m^ < 1, we find 

£® ~ \x{eBf 



1 + W, 



(B.5) 



Let us now consider the thermal contribution. Similar generalizations as in the vacuum 
case leads for // = to 



•(c) 2 



-"mat 



' ,; vv/./,. ' 



\2(7 



;b.6) 



Expanding the distribution function according to Eq. (|A.2|) , and using Eq. ( |A.3| ), we may 
again perform the Gaussian integral in p z , and sum the infinite geometrical series in n. 
Subtracting the contribution for B = 0, we find 



Ac) 

•'-'mat.l 



' 1 \2a oo 1 



1=1 



P 2 m}l 2 Jo t 2 



00 dt 



exp 



1 

2t 



\f3 2 l 2 eBt ,1 



smhaf3 2 l 2 eBt) 



exp(-(3 2 l 2 eB^st) - 1 



(B.7) 



We may again expand according to eBj3 2 l 2 t/2 < eB/rrip <C 1. Using Eq. ( A.4 ) the result 
reads 



C 



where we have defined 



(<0 



1 



uiai . j - ^x(eB) 2 [h{mp/T) + 0(eB/m 



I rT (x)=2^T(-l) 2 ° l xlK 1 (xl) 
1=1 



(B.8) 



(B.9) 



For small z, zK\{z) ~ 1, and for large 2, zK\(z) ~ y vr^/2exp(— 2) H43H , so that the sum 
is convergent. Now, to leading order, = m G oc T, so I a {mp/T) is T-independent. For 
completeness we may also write down the field independent contribution 



£mat,o(c) = m 4 



-1 



,2(7 



2tt 2 



EE 



\2cr/ 



^ /3m/ 



;b.io) 



We have found no analytical expression for this series, but its leading behavior is 0(m%T ), 
i.e. 0(g 2 T 4 ), sub-leading for small couplings. Of course, this field independent part is 
irrelevant for the effective coupling. 
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